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Abstract 

The mathematical problem concerning intrinsic storage optimisation is formulated and 
solved by means of variational analysis. The solution, though obtained in implicit form, still 
sheds light on many important features of the optimal exercise strategy. It is shown how 
the solution depends on different constraint types including carry cost and cycle constraint. 
Additionally, the relationship between intrinsic and stochastic solutions is investigated. In 
particular, we show that the optimal stochastic exercise decision is always close to the 
intrinsic one. 


1 Introduction 

The problem of storage optimisation is driven by the necessity of the storage owners to minimise 
their expenses and maximise the potential profit which can be gained by operating the storage. 
There are plenty of real world storage examples: Gas storage, Oil storage, Hydro power plant, 
Coal stock, etc. Another popular product on the gas and electricity markets is a swing contract, 
which is mathematically equivalent to the storage problem. By following a clever strategy - 
buying the underlying commodity cheap, storing it, and selling as the prices go up the storage 
owner can make a profit. Knowing when, and how much, to buy or sell determines the exercise 
rule for the storage. Thus, there is a need for exercise strategy optimisation - finding an exercise 
rule that maximises profit and minimises the risks. 

An exact storage option pricing method must also take the stochastic nature of the prices 
into account. This is usually done by modelling prices with a stochastic process that, as far as 
possible, captures the value-generating statistical properties of the real world price movements. 
This price model is normally calibrated on the available market data before it being used for 
the actual pricing. We refer to the corresponding option value as stochastic value. Alternatively 
one may consider the storage option in the case when the prices are fixed (to the forward prices 
observed on the valuation date). The corresponding option value is named “intrinsic”. The 
intrinsic value of the storage option is necessarily smaller than the stochastic value, since it does 
not take the time value into account. 

One of the most popular ways of storage optimisation is the “dynamic programming” al¬ 
gorithm. It can be used to optimize both the intrinsic and stochastic problems. Stochastic 
methods approximate stochastic prices by modelling either on a tree or by Monte Carlo simula¬ 
tion. A very efficient numerical algorithm for valuation using Monte-Carlo pricing was suggested 
by Longstaff and Schwartz 0]. A comprehensive description of the Least Squares Monte Carlo 
(LSMC) algorithm applied to the Storage Option is given in the paper by Alexander Boogert 
and Cyricl de Jong pQ. Dynamic programming offers the most generic approach, as it can deal 
with virtually any type of constraints. If the constraints are sufficiently simple then the intrinsic 
problem can be reformulated as a “linear programming” problem, which has much more efficient 
numerical algorithm, although usually the numerical efficiency is not an issue for the intrinsic 
problem. 

The numerical solution of the storage option provides answers to the two most important 
questions: what is the optimal exercise decision at every time moment given the current state 
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(prices and volume level of the storage), and what is the expected profit, provided we follow 
the optimal exercise strategy. Another important task is to find an optimal hedge strategy - 
a portfolio of derivatives (futures, options or any other financial instruments) - which would 
minimise the financial risks. Of course, market risk is only a feature of the stochastic problem, 
as the solution of the intrinsic problem is deterministic. Non-market risks, such as operating 
risks of a physical storage, are not taken into account as they represent systemic risk, and so are 
in practice determined ex post facto. 

In this paper we develop an analytical approach to the intrinsic storage optimisation problem, 
and consider a link between intrinsic and stochastic solutions. The solution is obtained in implicit 
form, and hence can not be used for the option value calculation directly. However it highlights 
many interesting features of the problem, and covers a large number of applications. One of the 
most remarkable results is related to the optimal exercise decision. We show that the intrinsic, as 
well as stochastic, exercise decisions are based on a so-called trigger price, and most importantly 
- the intrinsic and stochastic trigger prices coincide in the leading order. This implies that 
a correct stochastic exercise decision can be based on the intrinsic trigger price. Although a 
“rolling intrinsic” strategy is often used by traders, it is commonly believed to be sub-optimal. 
This work shows that the intrinsic exercise is in fact optimal, and we also demonstrate this with 
a numerical simulation. 

Our analysis also shows how different constraints can be modelled in an efficient way. Among 
the most interesting examples is a cycle constraint, which is sometimes imposed on storage 
contracts. Treating a cycle constraint numerically is very computationally expensive task, as 
it requires an additional dimension in the state space. We show that the cycle constraint can 
be mimicked by a virtual injection/withdrawal costs, which makes the numerical solution much 
simpler to implement, and is by orders of magnitude faster. 

A more detailed overview of the results is given in the section ’’Discussion and numerical 
analysis 11 . 

The paper is organised a follows. In Sec. [2] we give a formal mathematical formulation of 
the problem in terms of variational analysis. In Sec. [3] we present the solution to the intrinsic 
problem, and consider different special cases of constraints in Sec. [2 In Sec. [5] we briefly consider 
certain aspects of the stochastic problem, and give a proof that the optimal stochastic exercise 
is bang-bang, and can be approximated by the intrinsic exercise. A reader who wishes to skip 
the detailed mathematical derivations can go straight to the discussion of the results in Sec. [G] 

2 Problem formulation 

The storage problem can be formulated as follows: The storage option holder is given a right to 
store some amount of an underlying asset (for the current discussion let it be natural gas) in a 
(virtual) storage facility. At each moment in time the option holder may “do nothing”, inject 
gas into the storage, or release gas from the storage. Every time the gas is to be injected into 
the storage it must be bought on the market at the current spot price. Likewise, every time gas 
is released from the storage it is sold on the market. Since the market price of gas changes with 
time, this may lead to a non-trivial cash-flow 

The injection and release of gas may have to satisfy some operational constraints (for instance 
maximum injection and release rates, maximum storage capacity, etc.) specified as boundary 
conditions. Every exercise profile (trajectory in the time-volume space) yields a different profit. 
In this sense, the profit becomes a functional on the exercise trajectory, mapping each possible 
profile of injection and release decisions to a corresponding profit. The aim of the storage option 
holder is to maximise the profit by choosing an optimal exercise strategy. The problem can thus 
be formulated in terms of variational analysis - an optimal trajectory is the one delivering the 
maximum of the profit functional. 

If the market prices are static, that is, the prices for all delivery days are known in advance 
and never change, then the problem is deterministic. The possible profit is bound from above 
and from below, and thus there exists a trajectory such that no other trajectory yields higher 
profit. The maximal profit of the deterministic problem is called “intrinsic value”. We shall now 
investigate the deterministic problem by means of variational analysis. 
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We consider the storage problem in a continuous time approximation. Let T e be the total 
time of the storage contract, such that t £ [0,T e ]. 

Let F(t) be the market price of gas for delivery at the time t. The curve F(t) is called 
the forward price curve, since it is observed on the market at the time t = 0 and contains the 
information about the prices of gas with delivery in the future. By the definition of the static 
problem the forward curve never changes, and hence there is no difference between observation 
and delivery timcjj]. 

Let q(t) be the amount of gas contained in the storage at time t. We refer to the volume 
injected into (or released from) the storage, per unit time, as the “exercise”. The curve q(t) = 
dq/dt defines the exercise trajectory. The initial and terminal conditions are 

q( 0) = Qstart q(T e ) = Qend (2.1) 

We consider the case of the free terminal condition, where q(T e ) is not fixed, in Sec. 14.21 

The cash flow resulting from trading the gas according to the exercise strategy q[t) over the 
time interval dt is given by 


— q(t) F(t) dt 


( 2 . 2 ) 


An additional cash flow can occur from the costs of injecting or releasing gas (operating costs). 
Let us designate 7 (q) as the operating costs per unit time (we assume here that the operating 
cost is only a function of the exercise q). The terminal profit or loss is given by the cumulative 
cash flow over the storage lifetime. We thus introduce the target functional 


So[q(t)] 


q(t)F(t)+'y{q(t)) 


dt 


(?( 0 ) — Qstart 1 — Qend 


(2.3) 

(2.4) 


The value of this functional on the trajectory q(t) gives the storage value conditional on that 
trajectory occurring. Here we have used the subscript 0 to indicate the unmodified integral. In 
the next section we will introduce a modified integral S'[( 7 (f)], which includes additional terms 
intended to enforce the operational constraints. Looking ahead, we notice that on any trajectory 
q(t) allowed by the constraints the values of the modified and unmodified functionals coincide. 

Using “physics” terminology we introduce the (unmodified) Lagrangian Lq corresponding to 
the action functional Sq as 


Lo = ~q(t) F{t) - y(q(t)) (2.5) 

As mentioned above, the target functional is bounded both above and below, and hence there 
must exist trajectories on which the functional achieves its maximum and minimum values. We 
now make use of variational analysis to search for the extremal trajectory, i.e. the trajectory on 
which the first variation of the target functional vanished 

The types of operational constraints may differ depending on the different storage option 
types. Below we consider some typical constraints found in practice, which can be classified as 
local constraints. A local property implies that a constraint at time t can be expressed completely 
in terms of state variables and their derivatives q(t), q(t), F(t), F(t ),... at each time t. 

We consider the following two operational constraints: 

1 In reality the forward curve changes during the trading, and thus it is a function of both observation time t 
and delivery time T: F(t,T). 

2 We do not conduct a thorough analysis of existence and uniqueness of the solution of the variational problem. 
However, we point out that the class of functions q(t) on which the target functional is defined and achieves a 
maximum is rather broad. In particular, this class includes all continuous locally integrable functions. In some 
cases the solution space can even be extended to discontinuous or singular functions. We assume that all functions 
representing problem parameters - forward curve, boundary conditions, etc. - are sufficiently well behaved such 
that all mathematical expressions are well defined. If this is not the case, then we can apply the procedure of 
regularisation. We also utilise the concept of convergence of functions, which we will always understand to be 
weak convergence, i.e. ijjk —¥ i/j if S'fV’fc] —> *S’[y]. 
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( 2 . 6 ) 


1. The volume q(t) must remain within the interval 

Q{^) £ [Qminit')i Q max (*)] 

where the boundaries Qminit) , Qmaxif) can be time dependent. 

2. The injection/release rate q(t) is bounded by 

<i{t) e [ ^mini Tmax\ (2.7) 

Generally the maximal injection/release rates may depend on both the time and volume: 
r = r(t,q(t)). In the following analysis we will only consider constant injection/release 
rates. The case of volume dependent rates will be considered in Sec. 14.41 

One of the typical examples of a non-local constraint is the so called cycle constraint. It can 
be formulated as follows: One introduces an intake cycle variable as 

c (t) = j o iij) e{q{r)) dr where 9{x) = j * (2.8) 

which corresponds to the total injected volume up to time t. The cycle constraint requires that 
the terminal value c(T e ) does not exceed a certain threshold 

c{T e ) < Cmax (2.9) 

Similarly one can introduce a release cycle constraint, bounding the total released volume. We 
will consider the cycle constraint below in section f4.31 


2.1 Penalty functions 


To restrict q(t) from going beyond the range [ Qmin, Qmax\ , we can introduce a parametrised 
penalty function — 0[g(t), Nj,] and add it to the Lagrangian. The penalty function can be any 
smooth convex function, which in the limit —> oo approaches zero within [ QmimQmax ] and 
positive infinity otherwise. A particular example of this function could be 


a{q(t) - b) 


1 2 N+ 


where a 


2 v 


Qmax Qr 


b = 


Q 


max 


4" Qmin 
2 


1 



The parametrised representation has a purpose of regularisation , as the variational problem 
satisfies the strict requirements of classical variational analysis for finite N^. The transition to 
the limit —> oo should formally take place only after the finite N$ solution is found. 

A similar regularised penalty function —ip[q(t),N^,\ can be introduced to restrict q(t) from 
exiting the interval [r m i n ,r ma x\- In the limit —> oo it approaches zero within [r m i n ,r m ax} 

and positive infinity otherwise. The modified Lagrangian becomes 


L = L 0 - 0(g) - 0(g) 


qF + 0(g) + 0(g) + q(g) 


( 2 . 10 ) 


In the next section we find a formal solution to the variational problem, and show that the 
solutions of modified and unmodified constrained problems coincide. 


3 Solution of deterministic problem 

3.1 Euler-Lagrange equation 

We know from variational analysis that the extremal trajectory of the integral f L(q, q) dt satisfies 
the Euler-Lagrange equation (see e.g. 0 ) 

dL(q,q) _ d dL(q,q) 

dq dt dq 1 ' ’ 
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which is valid for both free and fixed terminal conditions. For the modified Lagrangian we obtain 
the following Euler-Lagrange equation 


<l>'(q) = Yt F + + 


(3.2) 


There are two types of solution for the latter equation. One is obtained on the interval where 
the trajectory remains strictly within the boundaries 

q(t) e {Q min i Qmax ) 

The second solution type is where the trajectory lies on the boundary 

q(t) = Qmin(t) Or = Qrnaxit) 


The general solution consists of pieces of both the solutions types. Let us consider each of these 
types separately. 


3.2 Solution between the boundaries 

First we consider the solution in the interval where the trajectory does not touch the boundary. 
In the limit N$ —> oo the l.h.s. of Eq. (ET21) vanishes. Thus we have 

F(t)+'ip , (q) + j'(q) =C (3.3) 


where C is constant for the whole period of time, where the solution does not touch the boundary. 

Let us consider one particular example of the operating cost function. It is typical for a 
gas storage facility that the operating cost is proportional to the amount of the gas released or 
injected into the storage. Again we can parametrise it with N . 7 so that for finite _/V 7 the function 
"f(q) is smooth, and in the limit iV 7 —> oo it is piecewise linear 


7 ( 9 ) 


7 inj q for q > 0 ; 

-7 rei q for q < 0 ; 


(3.4) 


where 7 i n j > 0 and el ^ 0 - 

The solution of Eq. (liOl) is now straightforward to find graphically (see Fig. [[]). In case 
Train < 0, r max > 0 we obtain: 


q(t) = 


T min 

0 


w 

=(*) 


F(t) > C + 7 re l 

C 'Yinj ^ A'(t) < C T 7 rel 
F(t) <C - 'Yinj 


(3.5) 


We see that the solution has three regimes, and the value C can be interpreted as a trigger 
price. If the price F(t) is within the “dead zone” [C — 'yinj, C + 7 re i\, the optimal exercise 
is “do nothing”: q = 0. If the price is below the dead zone, the gas is injected at maximal 
rate: q(t) = r ma x(t). If the price is above the dead zone, the gas is released at maximum 
rate: q(t) = r m i n {t). This result is a generalisation of the well-known Pontryagin’s maximum 
principle. 

In particular, without operating costs the solution has only two regimes injection and 
release (this feature is usually called in literature “bang-bang”). The transition between regimes 
occurs every time the forward curve crosses the trigger level C. 

The constant C has to be chosen in such a way that the solution q{t) satisfies the initial and 
terminal conditions 


^(O) — Q start 


q(T e ) — Qend 


Q 


start 



and thus can not be represented analytically in the general case. 


(3.6) 
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iP'(q)+Y(q) 


C — F(t) 


iinj 


+ 


T min 


OVeZ 


H-► 

T max Q 


Figure 1: Graphical illustration of Eq. ( 13 . 31 ) . 


3.3 Solution on the boundary 

There is not much we need to say about this part of the solution. On the boundary we obviously 
have 


q(t) = Qmax{t ) on the upper boundary 

q(t) = Qminit ) on the lower boundary 


There are two different points of view on the trajectory lying on the boundary. We may consider 
our problem as a general variational problem with a spatial boundary condition. A solution, 
that maximises the functional, can consist of parts going along the boundary and extremals - 
pieces of trajectory lying between boundaries and satisfying the extremal condition (i.e. having 
vanishing first variation). The Euler-Lagrange equation does not make sense on the boundary, 
since the trajectory can not be varied there. 

On the other hand, if we consider a regularised problem, the Euler-Lagrange equation does 
make sense on the boundary. We can derive the special case of the equation on the boundary. 
If the boundary Q m ax{t) (or Qminit)) is a sufficiently slow function of time, then the function 
ip'iq) vanishes, and the function <j>'{q) takes a finite value (positive if q = Q max or negative if 
q = Qmin)- Eq. ( 13 . 21 ) on the boundary becomes 




F{t) + 7 \q) 


From this follows 


Fit) + "f'iQmax) 
Fit) + j'iQmin) 


In particular for Q ma x = const and Q m i n 

Fit) > 0 

Fit) < 0 


>0 on the upper boundary; 
<0 on the lower boundary; 

= const we obtain 

on the upper boundary; 
on the lower boundary; 


(3.7) 

(3.8) 

(3.9) 

(3.10) 

(3.11) 


3.4 Approaching the boundary 

An important statement can be made about the connection point between extremals and the 
boundary. Let us first consider a simplified problem (a generalisation is straightforward), for 
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which we suppose a constant boundary condition, and zero operating cost: 

Qmax = const Qmin = COUSt 'j(q) = 0 

Let q{t) be the optimal trajectory and t* be the time when the trajectory touches the boundary, 
such that, for t < t* the trajectory satisfies the Euler-Lagrange equation and for t > t* the 
trajectory follows the boundary, e.g. the upper boundary: q(t > t*) = Qmax- 

Let us now consider a supplementary problem for the time horizon t £ [0, t*] with a terminal 
condition q(t*) = Qmax ■ We designate S sup the value of the supplementary target function. 
Obviously, the optimal trajectory of the supplementary problem coincides on [0, £*] with the 
original optimal trajectory 

Qsup(t)=q(t) t£[0,t*} 

Now consider a small variation of the boundary touch point 

t* = t* + St* 

For the supplementary problem we will obtain a new optimal solution q sup (t). As will be shown 
in Sec 13.6.21 the change of the value of the supplementary target function can be obtained from 
the time derivative of the boundary problem at the point t = t* , which is given by 

= Qsup (c - F(t*)) (3.12) 

SS sup = ^^St* (3.13) 

where q sup = q sup (t* — 0 ) = r max is the optimal exercise of the part of the trajectory just before 
the boundary touch. 

Regarding the original problem, we can now ask how the value of the target function has 
changed on the new trajectory q(t), which consists of q(t) for t < t* and is unmodified for 
t > ma x(t* ,t*) 

( q S u P (t ), t <t* 

q(t) = < q(t) t > ma x(t*,t*) 

{ Qmax t* <t< t* if t* < t* 

For the original variational problem the point t* is an internal point, and thus the new trajec¬ 
tory is a small variation of the optimal trajectory. Thus the original target functional remains 
unmodified: 

SS = 0 

The change of the value of the original target function is only due to the change in the extremal 
part of the trajectory for t < t*, since the part of the trajectory on the boundary does not 
contribute to the functional value. We thus conclude that the variation of the supplementary 
target function should vanish 

SS sup = 0 

Comparing this with Eqs. (13.121) and (13.131) we can conclude that, at the time of boundary touch, 

F(f) = C (3.14) 

In the general case (which includes time dependent boundary conditions and operating costs) 
we obtain the following conditions for the boundary touch points 

1. The trajectory touches the lower boundary from the left (Fig. [2ja): 

F(t*) = C + lrel , F(t*)< 0 (3.15) 

2. The trajectory touches the upper boundary from the left (Fig.[2]b): 

F(t*)=C- linj , F{t*)> 0 (3.16) 
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3. The trajectory touches the lower boundary from the right (Fig. [2]c): 


F(t*) = C - linj , F(t*)< 0 (3.17) 

4. The trajectory touches the upper boundary from the right (Fig. E]d): 

F(t*) = C + j rel , F(t*)> 0 (3.18) 



Figure 2: Boundary touch cases: a) lower boundary from left; b) upper boundary from left; c) lower 
boundary from right; d) upper boundary from right. 


These conditions must be satisfied for the time interval t £ (0,T e ). The extremal trajectory at 
the times t = 0 and t = T e does not have to satisfy these conditions as they are not internal 
points of the trajectory. 


3.5 Conclusion 

Eq. (13.511 together with the boundary conditions (13. 15113.lHll and a condition on the start and 
end volume EH provides an implicit solution for the problem. We summarise some properties 
of the solution: 


1. The optimal exercise strategy between boundaries is bang-bang, i.e. there are only 3 
possible optimal actions at each time moment: maximum injection, maximum release or 
“do nothing”. 

2. Generally the optimal trajectory consists of a number of sections - those following the 
boundary and interior sections satisfying Euler-Lagrange equation. 

3. For each piece of trajectory, separated from others by a boundary touch, there is a trigger 
price C and a “dead zone” [C — 7 i n j , C + 7 re i]- The volume in the storage is increasing at 
maximum rate if the price is below the zone, decreasing at maximum rate if the price is 
above the zone, and is constant if the price is in the zone. The width of the zone is defined 
by the operating costs. Parts of the trajectory separated by a boundary touch can have 
different trigger levels. 

4. If the trajectory touches the boundary in the time interval t £ (0,T e ), it must satisfy one 
of the boundary conditions (13.15H3. lHl) . 

5. Generally a trajectory satisfying all the previous conditions is not unique. Any solution 
must obey all these conditions, but not any function obeying all conditions is an optimal 
solution. 


The solution found in the previous section maximises the modified target functional S'[g(t)]. 
Let q(t) be the obtained solution. The optimal path q(t) can be substituted to the unmodified 
functional So [q(t )], on which we expect to obtain the value of the option. There are two questions 
we need to address: 1 ) if the optimal trajectory q(t) respects the constraints, and 2 ) if the 
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trajectory q(t), which is optimal for the modified functional, is also optimal for the unmodified 
one. For both these questions the answer is affirmative. In Appendix Owe show why. 

From this point onwards we drop the subscript “0” for the action integral, unless we want to 
emphasise the difference between modified and unmodified target functions. We also drop the 
“bar” for the optimal trajectory q(t) everywhere where it does not lead to ambiguity. 


3.6 Dependency on initial and terminal boundary conditions 

The optimal path q[t) and the target functional S'[q(t)] are calculated for fixed terminal point 
{t = T e ,q = Qend}- If this point is slightly shifted in either direction time or volume - the 
optimal path becomes different, and the target function as well. In this sense the target function 
can be viewed as a function of the terminal point S = S(q,t). In this section we find the 
derivatives of this function. 

From the standard variational analysis we have: 

dS_ _ dL 
dq dq 
dS _ .dL 
~dt ~ ~ q ~d^ 

We apply these equations to the modified objective function 5'[q , (t)]. The obtained formulas will 
also be valid for the unmodified objective function Sq [q(t)]. 


(3.19) 

(3.20) 


3.6.1 Volume derivative 


Using the definition (12.101) of the modified Lagrangian and Eq. (13.191) we find 


as 

dq 


F + i/j'(q) + 7 \q) 


The spacial derivative makes sense only within the boundaries. Making use of Eq. (13731) and 
replacing S with Sq we finally obtain 


dS 0 

dqend 


(3.21) 


where C is the trigger level of the last part of the trajectory. 

The derivative with respect to the initial condition q s tart is given by inverting the sign, 


dS 0 

dqstart 


(3.22) 


where C is the trigger level of the initial part of the trajectory. 


3.6.2 Time derivative 


From Eqs. (12.101) and (13.201) we find 


dS 

~dt 


q 

q 


i>'{q) + i{q) 
^\q)+i\q) 


</’(<?) - ip(q) - 7 (?) = 
7 (?) 


(3.23) 


Within the boundaries we make use of the relation 

9So 


dt en d 


= q 


. Replacing S with Sq we obtain 
C - F(t end ) 1 - 7 (g) (3.24) 


where q = q(t en d — 0) is the optimal exercise on the final part of the trajectory, obtained 
for a fixed terminal condition. In particular, one can see that dSo/dt en d > 0 (provided that 

Tmin ^ 0 T^max ^ 0 )- 
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On the boundary in the special case Q m ax = Qmin = 0 one obtains <7 = 0 , 7 (q) = 0 and 


dS 0 

d^end 


= 0 


(3.25) 


It is worth noting that the derivative with respect to the terminal (or initial) boundary condition 
is universal: it is independent on whether the optimal trajectory touches the boundary or not. 
To calculate the derivatives, one just needs to know the solution q(t) and the trigger price C at 
the end (or beginning) of the trajectory. This also remains valid with imposed cycle constraint 
ISec. 14.31) . 


4 Other constraint types 

So far we have only considered the optimisation problem with the simplest configuration. Often, 
realistic storage and swing contracts have additional constraints imposed, driven by operational 
or financial reasons. In this section we consider some typical additional constraints for storage 
and swing contracts, and discuss what influence they have on the solution. 


4.1 Carry cost 

Some storage contracts may include a “carry cost”, which is the cost of keeping the underlying 
commodity in storage. For instance, an oil storage is subject to a heating cost (fuel oil in a 
storage has to be kept warm at approximately 60°C). 

The carry cost is specified as a time-dependent price y c (t) per unit time per unit volume. 
Thus, the (unmodified) target functional becomes 


So — 


T e r 


Q p (t) +7(9) + 7e 


dt 


(4.1) 


Substituting the new Lagrangian into the Lagrange-Euler equation we obtain between bound¬ 
aries, 


p (t) +ip'(q) + i{q) = C(t) 


(4.2) 


where 


dC(t) 

dt 


7 c{t) 


(4.3) 


Thus, the solution of the problem with the carry cost is similar to that without, but has one 
difference: the trigger price is not a constant, rather, it is a growing function of time, satisfying 
Eq. (14.31) . 

It is easy to show that the carry cost can be replicated by an additional discount factor. 

Indeed, the trigger times occur when the forward curve crosses the trigger level F(t) = C(t). 

Using the solution for Eq. (14.31) in the form C(t) = Cq + f 7 c(w) du we find the condition for the 
trigger time 

F(t) - [ 7 c {u) du = Co 
Jo 

Thus the solution for the problem with carry cost is equivalent to the solution with modified 
forward curve F(t) = D(t) F(t ), where the discount factor D(t) is given by 

D (t) = 1 - Iciu) du « exp J q 7 c(u) duj (4.4) 


10 








4.2 Solution with free terminal condition 


Often the storage problem can be defined with a free terminal condition. For instance, one may 
be allowed to leave an arbitrary amount of the underlying in the storage, and for any remaining 
volume one gets a terminal pay-off, equivalent to selling that volume at an effective price F e , 
referred to as the residual unit price. We define the (unmodified) target functional S as 


rTe r 


S[q(t)] = - 


q(t) F (t) + 7(5) 


dt + Qend F e 


^(O) — Q start 


(4.5) 


where q en d — q(T e ). The optimisation problem is now to find an optimal trajectory q(t), satisfy¬ 
ing operational constraints, that would maximise the target function (14.51) . It is easy to see that 
any extremal trajectory of (ESI) is also extremal of (12.51 12.41) . Indeed, if the trajectory q(t) is an 
extremal of EHD , the variation of the target functional vanishes on any allowed small variation 
Sq(t). Since there is no terminal boundary condition, the variation Sq(t) may not vanish at the 
end point. But if the variation of the functional vanishes on the whole class of allowed variations 
Sq(t), it also vanishes on the sub-class of variations Sq(t) preserving terminal level ( 6q(T e ) = 0). 

We conclude that, in order to find a maximising trajectory for the functional ESP , we 
need first to find a solution for the functional (12.5112.41) with fixed terminal condition. Then, 
considering this solution as a function of terminal state q en d, we need to find a maximum of (14.51) 
as a plain function of q en d , provided q en d is not lying on the boundary. Thus, if the trajectory 
terminal point does not lie on the boundary, the target functional must satisfy the following 
condition 


dS 

dtqend 


= 0 


(4.6) 


This is valid for any target functional, regardless of whether the residual price is applied or not. 

As shown in Sec. 15.61 if the target functional El is considered as a function of the terminal 
state q en d, then its derivative with respect to q en d , in any interior point, is given by 


dqend 

Thus, the target functional (14.51) . considered as a function of q en d , has a derivative 


dS 

dqend 


= F e -C 


(4.7) 


where C is the trigger price of the final part of trajectory. We conclude that the condition for a 
target functional, with a free terminal state, to have a maximum is 


C = F e 


(4.8) 


which is valid if the terminal point is not lying on the boundary. 

4.2.1 terminal state cases 

We conclude that the terminal point of the optimal trajectory can only lie between the boundaries 
if C = 0 (assuming no residual unit price). If C > 0 then the trajectory inevitably terminates 
on the lower boundary. Similarly, if C <0, the terminal point lies on the upper boundary. The 
following relations must hold 


Qend — Qmin 


c> 0 

(4.9) 

^ Qend < Qmax 


<7 = 0 

(4.10) 

Qend — Qmax 

=> 

c < 0 

(4.11) 
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The inverse is also true with small modifications 


C > 0 

=> 

Qend — Qmin 

(4.12) 

C = 0 

=> 

Qmin — Qend — Qmax 

(4.13) 

C < 0 


Qend — Q max 

(4.14) 


If the residual unit price is paid, in the expressions (14.9114. 14|l we should substitute C —> C — F e . 

If the storage option has no residual unit price, it is most likely to terminate at the lowest 
possible volume level. Indeed, the upper boundary of the dead zone cannot be lower than the 
smallest forward price: 

C + el > min(T(t)) 

If the gas prices are positive, and operating costs are sufficiently small, then the trigger price is 
positive too. Hence the derivative of the option value with respect to the terminal level is nega¬ 
tive, and thus the optimal trajectory must terminate at the lowest possible level. Consequently 
the storage option can be considered as a problem with fixed terminal state (although it may 
not be fixed according to the contract). 

An opposite example is a swing option. A typical swing option is a contract between a gas 
supplier and a trading company. The trader buys gas at some predefined strike price from the 
supplier and sells it on the market at the current market price. The spread between the market 
and the strike price becomes the effective gas price for the trader. The swing option allows the 
trader to take the gas from the supplier according to some flexible scheme (i.e. to decide, when 
to take and how much within certain constraints). Thus, the swing option can be formulated in 
terms of a storage option. However unlike the storage option, the effective price (the spread) in 
the case of a swing option can be both positive and negative. Depending on the actual prices 
and on the contract constraints, it may happen that the terminal state does not lie on either 
the lower or upper boundaries. In this case we deal with the problem with the free terminal 
condition. For this problem the trigger level has to be equal zero. 


4.3 Cycle constraint 

A cycle constraint is sometimes imposed by the storage option seller in order to artificially restrict 
the option flexibility and reduce the option value. The cycle constraint is defined in (12.91) . For 
demonstration purposes we consider the intake cycle constraint. It caps the maximum injected 
volume during the operation period t £ [0, T e \. Instead of considering it in the form of inequality 
we can proceed as follows. The solution is found in two steps. On the first step we solve the 
problem without the cycle constraint and calculate the cycle variable on the optimal trajectory. If 
the cycle variable is below the threshold c morc , the obtained solution satisfies the cycle constraint. 
However if the solution violates the cycle constraint, we can impose a “biting” cycle constraint^ 

c(T e ) = f q(t) 9(q(t)) dt = c max (4-15) 

J o 

This condition allows us to formulate our problem in terms of conditional variational extrema. 
It is solved by means of a Lagrange multiplier, namely, we introduce a modified Lagrangian 


L = — 


qF + (f>(q) + ip(q) + j(q) + A q(t) 9(q(t)) 


(4.16) 


where A is an independent variable (the Lagrange multiplier). The solution to the latter problem 
will contain the undefined coefficient A, which can be found from the additional equation (14.151) . 

Within the boundaries, the Lagrange-Euler equation for the constrained Lagrangian becomes 
(compare with Eq. (13.31) 1 

C = F(t) + if>'{q) + 7 '(q) + A (e(q) + qS(q)^j 

= F(t)+ip'(q)+j'(q) + \9(q) (4.17) 

3 We say that the constraint is biting , when an inequality “<” becomes an equality. 
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Solving this equation we find that the solution is similar to the solution without the cycle 
constraint, but has a “dead zone” which is wider by an amount A. Thus, the cycle constraint is 
equivalent to having additional operating costs. One can find Active additional release and/or 
injection costs (constant for the whole trajectory), such that, if added to the original storage, 
will result in an optimal trajectory respecting the cycle constraint. It is inconsequential if the 
additional cost is for injection, release, or both, since it is only the width of the dead zone which 
matters. The only exception to this rule is if the optimal trajectory does not inject or release 
enough volume. For instance, if the trajectory never releases gas then the additional release 
cost will not effect the trajectory, and thus will not help respect the intake cycle constraint. 
Practically, one could use additional injection costs to fulfil the intake cycle constraint, and 
release costs to fulfil the release cycle constraint. 

Generally the solution between boundary touches contains two undetermined parameters - 
the trigger price C and the cycle Lagrange multiplier A. They can equivalently be expressed by 
the upper and lower bounds of the dead zone. These parameters have to be found in such a 
way that the obtained solution satisfies the boundary touch conditions (13.15113.181) . initial and 
terminal conditions (13.61) . and the cycle constraint (14.151) . 

Notice that the additional operating cost is the same for the entire time horizon, regardless 
of whether the trajectory touches the boundary or not. 

This result provides an efficient way of numerical computation of a storage problem with a 
cycle constraint. First one should solve the problem without any modification, and check the 
cycle variable. If the original solution breaches the cycle constraint, one needs to find a proper 
virtual operating costs, which allow the cycle constraint to be met. The search for the right 
virtual operating cost can be done by a “shooting” method, and in principle, should converge 
after a few iterations, since the cycle variable is obviously is monotonically decreasing function 
of the operating costs. 


4.4 Volume dependent injection/release rates 

If the maximal injection/release rates are volume dependent, the penalty function ip becomes an 
explicit function of volume: ip = ip(q,q). In this case the Euler-Lagrange equation within the 
boundaries becomes 


d q ip{q,q) 


d 

dt 


F + dqip(q,q) +j'{q) 


(4.18) 


The l.h.s. of this equation is an integrable implicit function of time. Designating 


C(t) = F + dq ip(q , q) + 7 '(q) 


we get 

=d q ip{q{T), q(r)); C(t) = C 0 + J d q ip(q(T),q(T)) dr 

We thus obtain the same solution as in the case of volume independent rates with the difference 
that C is not constant but a function of time. We conclude that the optimal exercise strategy 
still preserves the “bang-bang” property. However, it does not possess a constant trigger level. 
The function C(t) can be reverse-engineered from the optimal trajectory q(t), provided the latter 
is known. 


5 Optimal stochastic exercise 

In this paper we do not consider the stochastic problem in detail. However, we do wish to address 
one particular aspect of the storage option problem in its stochastic formulation, namely, the 
optimal stochastic exercise. 

In the stochastic formulation the prices are considered to be randomly changing during 
trading, as opposed to being fixed in the intrinsic problem formulation. The forward curve 
now becomes a function of two parameters - observation time t and delivery time T: F(t,T). 
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For every fixed delivery time T the forward curve becomes a stochastic process indexed by the 
observation time t. A choice of the price process poses a problem per se, and to some extent is 
driven by knowledge of the analytic solution properties. 

One peculiarity about the stochastic problem is that the optimal exercise profile can not be 
found upfront for the entire time horizon. It has to be continuously adjusted to the changing 
market prices. However for the current observation time t there will always exist an optimal 
stochastic exercise decision q s t(t), which remains valid for an infinitesimal time interval dt. The 
optimal exercise is conditional on the current state of the market F(t, T ) and of the storage q(t), 
or expressed in probabilistic language, conditional on the filtration Tt- We refer to the optimal 
exercise decision at time t as the optimal prompt exercise (the prompt exercise is a continuous 
time analogue of the spot trade in the discrete time systems). 

In a similar way we can introduce an intrinsic prompt exercise. For every time t we can 
formulate a new intrinsic problem, for which q(t ) plays a role of the initial condition, and 
F(t,T) is the forward curve (as a function of T). The solution of this problem is an optimal 
intrinsic exercise qi n {t, T), for which observation time t plays a role of a parameter, and the time 
derivative is taken with respect to T. The intrinsic prompt exercise is defined as the initial value 
of the optimal exercise profile qi n (t,t). Below we use the notation 

r{t) := q in (t,t) 


for the intrinsic prompt exercise. 

In this section we establish a link between the optimal stochastic q st ( t ) and intrinsic prompt 
exercise r(f), and show why the intrinsic exercise is a good approximation for the stochastic 
one. We do not consider injection/release costs for the sake of simplicity of notation. The 
generalisation is straightforward. 

Although the main result of the current section turns out to be independent of the choice 
of price process, for the sake of completeness we discuss some of most common approaches to 
modelling the stochastic price evolution. 


5.1 Forward curve evolution model 

Usually the forward price evolution is approximated by a zero mean stochastic process governed 
by the stochastic differential equation of the form 


dF(t,T) 

F(t,T) 


dM t (T) 


(5.1) 


where M t (T) is a martingale process of time t, parametrised by the variable T. For the derivation 
of the optimal stochastic exercise the particular choice of the price process is not essential. The 
only important assumption about the price process is that it is drift-less, i.e. the mean of the 
price increment (for every delivery time T) is vanishing: 


(dF(t, T)) = 0 


(5.2) 


This assumption is plausible on liquid markets without friction. 

For different delivery times T the processes M t (T) may or may not be correlated. One 
particular example of the price process is the 1-factor model 


dF{t, T) 
F(t,T) 


a dW t 


(5.3) 


where W is the standard Wiener process. The model was first presented in [5], its application 
to the Energy market can be found in [2]. The one factor model is very popular due to its 
simplicity. Notice that the processes M(T) for all delivery times T are perfectly correlated. 
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5.2 Stochastic trigger price 

Let F(t, T) be the forward curve observed at time t. We define the spot price process s(t) as 

s(t):=F(t,t) (5.4) 


Since the forward price curve evolves with time, the spot price s(t) = F(t,t) generally does not 
follow the original forward curve (s(t) ^ F(0,t) almost everywhere). The spot process s(t) is 
uniquely defined by the forward price process. 

Consider a set of possible spot price curves Si(t), each having corresponding probability pi. 
Viewing a particular spot price curve Si(t) as a deterministic function, we can solve a static 
optimisation problem for the target functional 


Si = 



q(r) Si ( t ) dr 


defined for the time moment t = 0. For every spot price curve Si(t) we can find a corresponding 
trigger price C\, optimal exercise strategy qi(t) and intrinsic value Si. The prompt intrinsic 
exercise 

n := qi( 0 ) 

for the time moment t = 0 is optimal on the i-th spot price path. 

Now we are searching for an optimal prompt exercise decision q st (0) for the stochastic problem 
at time t = 0. We omit the argument 0 in what follows to simplify the notation. The set of 
the spot price paths with corresponding probabilities contains the entirety of the stochastic 
information required for making an optimal stochastic exercise decision. Since we can only make 
one exercise decision at a time, this exercise decision will not be optimal for some spot price 
paths. If on the ith spot price path the optimal exercise is r l . the sub-optimal exercise decision 
qst 7 ^ i~i will lead to a loss of value on that path. We need to make such a choice of q st that 
would minimise the expected losses over all paths for which the decision is not optimal. 

During the time interval dt the volume increment in the storage is given by dq = q s t dt. Thus, 
on the ith spot price path the volume in the storage is by (q s t — ri)dt bigger than it would be if 
the exercise was optimal for that path. This leads to a change of the value on the ith path by 


5Si = (Ci - F(0,0)^ (q st - r^j dt (5.5) 

Indeed, the change of the value occurs for two reasons: the additional expense SSn of buying 
additional volume, and an increase in value SSi 2 due to a change in the actual storage volume. 
The price of additional volume is 


SSn = -F( 0, 0) (q st - r^j dt 

The change in value due to a change in the initial volume follows from the formula dS/dq s t. a rt = C 
(see Sec. 13.fill : 

SS i2 = Ci (q s t - r^j dt 

Combining the latter two equations we obtain Eq. (1531) . Now, averaging SS over all paths, we 
have 


(SS) = ^/Pi (Cj - F( 0 , 0 )) (q st - r^jdt (5.6) 

i 

4 The discussion of why the expressions for the target functional are well-defined would lead us too far from 
the scope of the paper. However we should point out that the problem can be formulated in a discrete time, 
where it becomes a well-defined finite-dimensional problem. The transition to the continuous time can be done 
if we demand that the spot price process is sufficiently “good”. Wiener process is known to be almost surely 
continuous, and the integral of Wiener process is well-defined. 
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where pi is the probability of the *th path. The optimal exercise decision q st should maximise 
the expected change of option value (SS). Deriving (SS) with respect to q s t we obtain 


^7 = E Pi ( Ci - F (°-°)) dt = ((C)-F( 0,0)) dt 


(5.7) 


where (C) = JA Pi Ci is the trigger price averaged over all possible spot price processes. If the 
spot price F(0, 0) is bigger than the average trigger price (C), then the derivative is negative, and 
the optimal stochastic exercise is the smallest possible allowed by the constraints. Similarly if 
F( 0,0) is smaller than the average trigger price, then the optimal exercise is the biggest possible 
allowed by the constraints. Thus, the optimal prompt stochastic exercise is 


qst = 


T min 
T max 


F(0, 0) > (C) 
F(0, 0) < (C) 


(5.8) 


We see that the optimal stochastic exercise is again bang-bang, and that the expected value of 
the intrinsic trigger price (C) can be interpreted as the stochastic trigger price C st : 


C st = (C) 


(5.9) 


Next we show that, under some assumptions, in the leading order of the Taylor expansion 
the stochastic trigger price equals the intrinsic trigger price. Indeed, let us designate so(t) as 
the spot price path which coincides with the initial forward curve at time t = 0: 


s 0 (t) = F(0, t) 


The trigger price Co is then the intrinsic trigger price, calculated at time t = 0 on the basis of 
the forward curve F(0,t). All other paths Si(t) are different, and we designate 

Ssi{t) = Si(t) - s 0 (t) 

This difference is zero at the beginning of time period: <5sj(0) = 0 since all paths of the spot 
price process start in the same point Sj(0) = F( 0, 0). At any future point in time the difference 
6si(t) may be arbitrary large, but the majority of paths stay within the standard deviation 

VWm 


If 5si(t ) remains small, then the difference in the trigger level SCi = Ci — Co can be calculated 
by means of the first order Taylor expansion (see App. 0: 


where we designated 


SCi 



S[Cq - s 0 (t)] Ssi(t) dt 


At* — T* max 


Tmin 


K = Ar I S[Co — so(t)] dt 

Jo 


Since the forward price process is a martingale, we conclude that (Ssi(t)) = 0, and hence in 
the first order (SCi) = 0. Consequently 


C st = (C) =C 0 + (SC) » Co 


(5.10) 


which means that the stochastic trigger price equals the intrinsic trigger price in the first order 
of the Taylor expansion. 

We performed the calculation of the stochastic trigger price for the initial time moment t = 0. 
For any later time we can repeat the same calculation, and thus the following statement holds: 
for any time t the optimal stochastic exercise is bang-bang, and the stochastic trigger price C s t(t) 
is approximately equal to the intrinsic trigger price C(t), calculated at time t on the basis of the 
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deterministic forward curve F(t,T) and initial condition q(t). The exact stochastic trigger price 
can be found as an average over the intrinsic trigger prices calculated for the entire set of spot 
price paths having F(t,T ) as initial condition. 

Note that to derive the relation (EU we did not use any approximation, and hence this 
relation is exact. In order to derive relation (15.1011 we used only the assumption of small variation 


So(t) 


< 1 


(5.11) 


For the standard 1-factor price process ([OP the relative spot price variation can be easily 
estimated, 


/exp 

^U-e- 2at ) 

- 1 

/ 

2a v / 



So(t) 

In particular, for aT e <l the condition for the intrinsic exercise approximation becomes 

cr 0 \fr~e < 1 


and for aT e >l we obtain the condition 


co 


a /2 a 


<C 1 


5.2.1 Equivalent stochastic trigger prices 

The intrinsic trigger price is a well-defined value, which has an impact on the whole intrinsic 
exercise strategy (or at least for the whole part of trajectory between boundary touches). One 
can interpret the trigger price as an alternative way of representing the exercise strategy, since 
the optimal exercise can be derived for any point in time from the trigger level and price. 

To the contrary, the stochastic trigger price only has an impact on the prompt exercise 
decision over the infinitesimal time interval dt. The trigger price may be different for every 
observation time, since, due to the stochastic nature of the prices, there is no optimal exercise 
strategy that can be calculated for the whole time horizon of the storage contract. 

The stochastic trigger price is thus not uniquely defined. Indeed, if some price C{t ) is a 
trigger price at time t, then it defines the prompt exercise according to the rule: inject if the 
spot price s(t) is below C(t), and release if the spot price is above C(t). In the first case, any 
price C(t) which is above the spot price would lead to the same exercise decision, and thus play 
a role of the trigger price. Similarly, in the second case, any price C(t) < s(t) would lead to 
the same prompt exercise. Thus we can say that for each time moment there exists a class of 
equivalent stochastic trigger prices. 

If C\{t) and C 2 (t) are two equivalent stochastic trigger prices, and s(t) is the spot price, then 
obviously 


s(i)<Ci(i) s{t)<C 2 {t) (5.12) 

s(t) > C\(t) <t=> s(t) > C 2 (t) (5.13) 

Form these statements we conclude that there is one case when all equivalent trigger prices 

should coincide 

s{t)=C 1 {t) ^ s(t) = C 2 (t) (5.14) 


5.2.2 Value loss due to sub-optimal exercise 

Next we show that, similar to intrinsic case, the stochastic trigger price (one from the range of 
equivalent prices) is related to the stochastic option value V as 


C st 


dV 

dq 
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where q is the volume level. 

Let V ( q , t) be the stochastic option value at time t, considered as a function of the volume in 
storage q and time t (for a given forward curve F(t,T)). Let the spot price be s = F(t,t). We 
consider a portfolio composed of the storage option (over the remaining time horizon), hedge 
trades and cumulative cash-flow generated by the spot exercise trades. The increment of the 
portfolio value cflL over the time interval dt is given by 

dV 

dU = —s dq + —— dq + ... (5.16) 

oq 

The first term is due to the cash flow from purchasing volume dq at the spot price s, and 
the second is due to changed option value as a function of volume q. We omitted the terms 
originating from derivatives with respect to time and forward price, as they have no impact on 
the derivation of the trigger price. The portfolio value increment is maximised by the following 
choice of volume increment 


f dqmax — Tmax dt s < dV/dq 

\ dqmin — ^min dt s > dV/dq 


(5.17) 


We see that the derivative dV/dq can be interpreted as a stochastic trigger price. It is equivalent 
to the trigger price C s t calculated previously. 

The trigger price allows us to estimate a loss of portfolio value due to a sub-optimal exercise. 
Considering the portfolio increment dll as a function of the exercise dq we can write 


ddU 

~dd<q 


= C st — s 


(5.18) 


Let dq be an optimal exercise, and dq be some arbitrary (sub-optimal) exercise, which leads to 
the portfolio increment dll. Designating Sq = dq — dq as the deviation of the prompt exercise 
from the optimal one, and <511 = dll — dll as the corresponding value loss, we obtain 


<511 = (C st (t) - s(t)) Sq 


(5.19) 


Notice that the portfolio value is insensitive to the sub-optimal exercise at the trigger times t* 
defined by the equation C st (t*) = s(t*). 


6 Discussion and numerical analysis 


In this paper we presented a solution to the intrinsic storage optimisation problem. It is repre¬ 
sented in an implicit form, as the solution contains an undefined parameter - the trigger price 
- which can not be expressed analytically. However, the solution gives us a deep insight into 
the essence of the storage problem, and provides evidence for sometimes very counter-intuitive 
features of the storage option. 

Here we first give an overview of the formal solution, and then discuss the most important 
conclusions that can be drawn from the analytic results. 


• An optimal exercise of the intrinsic problem can be expressed by means of the trigger price. 
This is such a price level C for which the optimal exercise decision can be represented as 


f r max (t) F(t) < C 

\ r min (t) F(t) > C 


( 6 . 1 ) 


If injection and release are subject to operating costs then there is also a “dead zone” 
around the trigger price. The optimal exercise in this case is represented in Eq. (13.511 . It 
contains an additional “do nothing” state, which is triggered when the price is within the 
dead zone. The width of the dead zone is solely defined by the operating costs. 

In the more general case, an optimal exercise trajectory may partially go along the upper 
or lower boundary. The entire trajectory then consists of several parts of “active” injec¬ 
tion/release, separated by a boundary touch. The trigger prices for different “active” parts 
separated by a boundary touch may be different. 
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• If the optimal exercise trajectory q(t) touches the storage boundary (i.e. reaches the storage 
minimum or maximum volume), then conditions (I3.15M3.1B1) must be satisfied. In the case 
of no operating costs these conditions are simplified to 

F{t*) = C (6.2) 


where t* is the time of the boundary touch. 


• If the option value V is considered as a function of the terminal level q en d , then is has a 
simple relation to the trigger price given by 


dV 

&Qend 


= F e -C 


(6.3) 


where F e is a residual unit price (price paid back for the gas remaining in the storage after 
the end of the contract). 


Now we are going to discuss some of the most significant practical implications, which follow 
from our analysis. 


• An important statements can be made about the terminal level of the storage. First 
consider a storage option with no residual unit price. If the prices are strictly positive 
and if operating costs are sufficiently small, then the trigger price is positive, and hence 
the derivative dV/dq en d is negative. Thus the optimal trajectory has to terminate at the 
lowest possible storage level 


Qend — Qmin (6.4) 

With a good accuracy the same applies to the stochastic storage option. 

• Unlike storage, for a swing option it is possible to terminate between global volume con¬ 
straints. Indeed, an effective price of the underlying of a swing contract is a spread between 
the market price and the strike price, which can be positive or negative. Hence the sign 
of the trigger price is indefinite. The actual terminal volume (in other words, the total 
taken volume) depends on the contract “moneyness”. If the swing contract is deep in the 
money (spread is positive), the trigger price is positive and the exercise trajectory will 
terminate on the lower boundary. If the contract is deep “out of the money”, the spread 
is negative, and so is the trigger price. Hence the exercise trajectory will end on the upper 
boundary. For an “at the money” option the total taken volume might lie between the 
global constraints. In this case the trigger price must vanish 


C = 0 (6.5) 

The same applies to the storage option, if a terminal unit price is paid for the remaining 
volume. Depending on the moneyness, the terminal volume may lie between the volume 
constraints. 

Notice that a vanishing trigger price implies a simple exercise strategy: release as soon 
as market price is above the strike. In this case the swing contract becomes a strip of 
independent vanilla call options. If the terminal point of the intrinsic exercise trajectory is 
sufficiently far from the global constraints, then the same applies to the stochastic swing 
option: it can be priced as a strip of independent vanilla call options. The independence 
follows from the fact that the majority of stochastic exercise trajectories will also terminate 
between the boundaries. Hence the trigger price C = 0 is preserved, and the exercise rule 
for all maturities remains the same. A swing option which can be represented as a strip 
of independent options has a significant advantage, since it can be replicated (and hence 
priced and hedged) by the market observables (of course, only to the extent allowed by the 
market liquidity). 
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• If a carry cost is imposed, then the trigger price becomes a growing function of time 
satisfying 



( 6 . 6 ) 


where 7 c (t) is the carry cost per unit volume per unit time. As has been shown in the 


main text, in terms of exercise strategy, carry cost is equivalent to an additional discount 
factor 



Notice that the discount factor leads to a correct exercise strategy, but wrong cash flow. 
The correct cash flow has to be calculated after the optimal exercise is found. 

• A cycle constraint is equivalent to an effective additional injection/release cost. Practically 
this means that in order to meet the cycle constraint one can use the following iterative ap¬ 
proach. First solve the intrinsic problem without any modification. If the cycle constraint 
is breached, one should add an operating cost (let us call it virtual opex) and solve the new 
modified problem. Too low a virtual opex might not be sufficient to meet the cycle con¬ 
straint, too high a virtual opex will reduce the cycle variable below the constraint. After 
a few iterations one should converge to an optimal choice of virtual opex, which leads to a 
biting cycle constraint. Such an iterative solution is computationally much more efficient 
than a brute force approach, since the latter would require an additional dimension in the 
state space. 

Notice that, similar to the case with carry cost, the virtual opex gives rise to a correct 
exercise strategy, but wrong cash flow. The correct cash flow has to be calculated after 
the optimal exercise is found. 

• One of the most important results of our analysis concerns the stochastic spot exercise. We 
have shown that the optimal stochastic exercise is “bang-bang” in exactly the same way 
as it is for the intrinsic problem. Thus the stochastic exercise can be expressed in terms 
of the stochastic trigger price C s t(t). Unlike the intrinsic problem though, the stochastic 
trigger price cannot be defined upfront for the entire time horizon. For each time moment 
t: it can only be determined for a short time interval, over which the forward curve does 
not change significantly as a function of the observation time. Typically gas storage and 
swing contracts are exercised with a daily granularity (day-ahead nominations), and hence 
one day is a natural time interval of the trigger price validity. A more interesting use case 
for the trigger price is posed by weekends and bank holidays, when there is no gas market 
activity, and the trigger price should be prolonged over a period of a few days. 

• The stochastic trigger price turns out to be very close to the intrinsic trigger price, or 
equivalently, the stochastic prompt exercise is almost identical to the intrinsic prompt 
exercise. The importance of this fact is difficult to overstate, as it justifies the so called 
rolling intrinsic (RI) strategy. This strategy is a rule to exercise the storage option using 
the intrinsic optimisation on every time step. The RI strategy is performed as follows. At 
each time step t , given the current market prices F(t, T) and the current storage level q(t), 
one solves the intrinsic problem and finds the optimal intrinsic exercise profile q(t,T) for 
the remaining time horizon. The prompt exercise (the exercise trade) q{t : t) is traded on 
the spot market, and the remainder (the hedge trade) on the forward market (of course, 
when performing exercise and hedge trades, one has to take into account the existing hedge 
position). This action is repeated each trading day until the end of the storage contract. 

The rolling intrinsic strategy is one of the most frequently used by traders. This strategy 
is attractive for several reasons, one of which is that it is trivial to follow. In particular, it 
does not require any knowledge of, or calibration to, the price process. 

Since the intrinsic hedge is not perfect (as it does not take the stochastic price dynamics into 
account), the terminal portfolio value is to some extent random, and may be significantly 
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different for different price evolution scenarios even with identical initial conditions. An 
important fact is that, on average, the terminal value of the portfolio is only determined by 
the prompt exercise^. Because the intrinsic prompt exercise, as we have seen, is very close 
to the correct stochastic prompt exercise, we conclude that the RI strategy on average 
generates the correct terminal option value (in a risk neutral probability measure). 

• A correct stochastic option value can only be achieved by following an optimal stochastic 
exercise strategy. However there are a few time moments during the storage contract life 
time, when the prompt exercise decision is irrelevant. Indeed, according to Eq. (15.191) the 
deviation of a portfolio value from its optimum is given by 

<HI = ( C st (t ) - s(t)) Sq (6.7) 

where Sq is the deviation from the optimal exercise, and s(t) is the spot price. We see that 
the portfolio value loss (5n becomes insensitive to the wrong exercise every time the spot 
price crosses the trigger level, i.e. at trigger times. 

An important consequence follows from this statement. Since the stochastic and intrinsic 
trigger prices are very close, they may only lead to a different exercise decision at those 
time moments when the spot price gets close to the trigger price, in other words - at 
trigger times. However, if instead of a correct stochastic exercise one performs an intrinsic 
exercise, it won’t lead to a change of option value, since the option value is insensitive 
to the exercise decision at the trigger times. This gives yet another justification why the 
intrinsic exercise is a good approximation of a correct stochastic exercise. 

Another interesting consequence is that, according to the boundary touch conditions (15.151 
Km . every time the trajectory is approaching the boundary, the option value becomes 
insensitive to the exercise decision. 

We performed a numerical simulation in order to demonstrate the effect of the intrinsic ex¬ 
ercise on the PnL. We simulated a PnL evolution of a portfolio consisting of a 1-year storage 
option, a linear hedge and a cumulative cash flow. For the simulation we generated 300 dif¬ 
ferent evolution scenarios of a forward curve, driven by the one-factor price process m with 
parameters a = 20 and cr = 0.9. 

For each price evolution scenario we ran two instances of the portfolio, one of them using a 
correct stochastic exercise decision as a benchmark, another using an exercise decision based on 
intrinsic valuation at every time step. 

Since the storage pricing model used the same assumptions as those used for simulating a 
forward curve, the first portfolio is perfectly hedged, and hence its PnL follows a straight line 
plus some noise which is due to the discrete time and Monte Carlo valuation error. 

The second portfolio used correct stochastic hedging, but intrinsic prompt exercise. The PnL 
of the second portfolio may be slightly different only due to the rare events when intrinsic and 
stochastic exercise decisions are different. Our aim was to demonstrate that the PnL of the two 
portfolios is essentially the same. 


Full option value 

305.8 ±0.5$ 

Time value 

31.1 ±0.5$ 

Terminal PnL (stochastic exercise) 

306 ±0.2$ 

Terminal PnL (intrinsic exercise) 

305.4 ±0.2$ 

PnL loss (intrinsic exercise) 

0.6 ± 0.1$ (2% of the time value) 


Table 1: Simulation results: intrinsic trigger price yields correct exercise decision. 

The results of the simulation are presented in the table [1] (we used the standard error for the 
accuracy interval). We found that the intrinsic exercise led to a systematic PnL loss of about 
2% of the option time value. 

5 Indeed, the hedge trades consist of pure financial instruments. Since the price process is a martingale, the 
average value of all hedge contracts vanishes. It is only the “physical” prompt exercise trades which generate a 
non-trivial average cash flow. 
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For our example we deliberately chose such parameters (fast storage and fast oscillating 
forward curve) to enforce a large number of trigger times, in order to make the effect of intrinsic 
exercise more pronounced. More realistic examples of the storage option would have much fewer 
trigger times, and hence even smaller PnL loss. Thus the simulation confirmed that the intrinsic 
exercise has a negligible effect on the PnL. 
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Appendices 

A Modified and unmodified functionals 

Let q(t) be the optimal solution of the modified functional 5[g(t)]. Here we show why q(t) 
respects the constraints of the original problem, and why q(t) is also an optimal solution for the 
unmodified functional S'o [</(£)] with constraints. 

The optimal trajectory q belongs to the broad class of locally integrable functions C\. Let 
us introduce a class of trajectories 1C C L\ satisfying all the constraints. By definition, penalty 
functions vanish on 1C: 

4>{q) = 0 ip(q) = 0 iff q £ K 

Obviously the values of modified and unmodified functionals coincide on 1C: 

%] = « V 9 £/C 

Next we show that the optimal trajectory q belongs to the class 1C and hence satisfies the 
constraints. Indeed, generally on some parts of the extremal trajectory ip'(q) 7 ^ 0 and (p'(q) ^ 0. 
According to Eas. (13.31) and (13.71) they take finite values. One can show (we leave it without 
proof) that in the limit —> 00 everywhere where <p' (or ip') take finite values, the function 

<p (or ip) vanishes. We conclude that on the extremal trajectory 

cp(q) = 0 ip(q) = 0 (A.l) 

and hence q £ 1C. Form (IA.1I) follows that the values of modified and unmodified target func¬ 
tionals coincide on the maximal trajectory: 

S[q\ = S 0 [ 9 ] 

The last step is to prove that the trajectory q delivers the maximum for the unmodified 
functional on the space 1C. It is obvious from the following consideration. The trajectory q 
maximises the modified functional on the space L\ and hence it will also maximise this functional 
on the smaller space 1C C L\. Since the values of both functionals coincide on /C, the trajectory 
q will also maximise the unmodified functional So on the space of all trajectories satisfying the 
constraints. 


B Variation of the trigger price 


Here we calculate how the intrinsic trigger price is modified by a small variation of the forward 
curve. We only estimate the first order Taylor expansion in this section. 

Let SF(t) be the variation of the forward curve, and SC is the corresponding variation of the 
trigger price. In order to derive the relation between SF(t) and SC, we depart from the relation 
connecting the initial and terminal conditions: 


Qend — Qstart T / 9’(^) dt 


(B.l) 
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where q(t) is the optimal intrinsic exercise strategy. It can be represented in the form 


q(t) = r m i n + A rO[C - F(t)] 


(B.2) 


where 


Ar — T max 


T min 



for x > 0 
for x < 0 


We simplified the expression for the optimal exercise by excluding the boundary touch events. 

Since the terminal level has to be preserved, the variation of both sides of Eq. El must 
vanish. Substituting Eq. (IB. 21) into Eq. (IB. II) and finding its first variation with respect to C 
and F, we get 


Ar (<f[C7 - F(t)] SC - 5[C - F(t )] SF{t)\ dt & 0 (B.3) 

where <5[a?] = 0'[x\ is the Dirac delta function. From here follows 

A /*X' e /*Te 

SC ~~FF S[C - F(t)} SF(t) dt where K = Ar 6[C - F(t)] dt (B.4) 

F J o Jo 
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